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INTRODUCTION 


Recently a new low-cost method has been investigated as a means of pre- 
paring very pure silicon from silane (SiH^K The fundamental chemical reac- 
tion of this technology is 

SiH 4 (gas) heat 9- Si I + 2H 2 (gas) . 

When cold silane gas is heated, it decomposes forming solid silicon* 

A number of methods of effecting this decomposition reaction are shown in 
Figures 1(a) and 1(b)* Figure 1(a) is a free space reactor in which cold 
silane is injected into an externally heated reactor* With this system there 
are problems associated with the consolidation of the extremely fine dust 
which is formed, and with the formation of large agglomeration or scabs of 
solids which grow on the heating surfaces* Another way of preparing very pure 
silicon is to fluidize a bed of hot silicon particles with cold silane gas * 
Figure 1(b) is a normal externally heated fluidized bed with a silane feed* 
Because the silane will tend to condense out on the hottest surfaces, the 
problem with this system is to keep the bed particles hot while keeping the 
walls of the reactor cool* Figure 1(c) is an internally heated bed wherein 
the heater is located far above the inlet point at which the SiH^ is expected 
to decompose* Circulation of the fluid bed particles is designed to keep a 
sufficiently high temperature of particles near the inlet so as to obtain 
complete reaction at the bottom of the reactor* 

The purpose of our research is to explore a still different processing 
method which uses radiant heating of the particles of a fluidized bed, as 
shown in Figure 1(d) * The process involves heating the particles at the 



2 


Free space reactor 


Fluidized bed reactor Internal heated 
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Radiantly heated 
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(b) (c) (d) 


Fig. 1 Different ways to make very pure silicon from silane. 
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surface of a shallow fluidized bed by radiant heaters located above the sur- 
face of the bed. The radiant sources could be cooled by a purge of inert gas 

which is exhausted with the spent The hot bed particles would circulate 

down to the bottom of the bed where they would heat the cold silane gas rising 
through the distributor plate. Decomposition and deposition of the silane 
then follow. 

Two important systems characteristics need to be evaluated. One is the 
effective absorptivity of the bed. This is the fraction of electrical energy 
supplied to the heaters which is actually absorbed by the particles . This 
involves the geometry and directionality of the heaters , the absorptivity of 
the surface of the fluid bed, and the flow of energy from the top to the 
bottom of the bed. This last named characteristic is primarily dependent on 
the circulation rate of the solids and the bed configuration, and is known to 
be very rapid. 

The second important factor is the heat transfer coefficient between the 
hot bed and the cool distributor plate. Knowledge of this will allow predic- 
tion of the temperature of the distributor plate during silicon production. 

It is important that the distributor plate be kept below about 350°C [1]; 
otherwise silane will decompose while passing through it , depositing silicon, 
and plugging it. Thus knowledge of the heat transfer coefficient will deter- 
mine what conditions (vessel geometry, feed gas temperature) will give satis- 
factory designs . 
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EXPERIMENTAL 


Apparatus and Materials 

An outline of the experimental equipment is illustrated in Figure 2. 
Figure 2(a) shows that the bed may be considered to be composed of three 
parts* The upper expanding section is a square section tapered from 394 mm to 
152 mm in a vertical distance of 440 mm* The lower section is a square sec- 
tion tapered from 152 mm to 13 mm with a taper angle of 20° . The middle 
section is a square section with a 152 mm square distributor plate* Figure 
2(b) is a square tapered bed with a 76 mm square distributor plate* Figure 
2(c) is a square tapered bed without distributor plate— we call this the 
square pyramid bed* Figure 2(d) is a normal 152 tom circular bed® It is made 
of 152 mm pyrex tubing with a 1000 W electric resistance heating coil just 
above the surface of the bed* 

For the “cold" tests the heater used was a 1000 watt , 120 volt , spotlight 
(Sylvania Q1000PAR64NSP) mounted in a water-cooled copper shroud* For “hot" 
tests, four such lamps were mounted in a water-cooled heater assembly* 

The distributor plate was made of a 1*6 mm thick sintered stainless steel 
plate with nominal 5 ym pores obtained from Mott Metallurgical Corporation* 

Silicon particles for the fluidized bed were obtained from Union Carbide 
Corporation and were measured to have the following size distribution: 


Mesh 


(U. S • Standard) 

Opening (ym) 

^screen, i 

(Wt%) 

+40 

425 

425.0 

0.02 

-40 +70 

425/212 

318.5 

33.04 

-70 +80 

212/180 

196.0 

27.35 

-80 +100 

180/150 

165.0 

28.54 

-100 +140 

150/106 

128.0 

10.24 

-140 +200 

106/75 

90.5 

0.55 

-200 

75 

75.0 

0.25 




5 
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tapered bed 
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Square Round bed 

pyramid bed 
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Fig. 2 Different geometries of the experimental beds used in this study, 

(a) 152 mm square bed with distributor plate * 

(b) square tapered bed with 76 mm square distributor plate, 

(c) square pyramid bed without distributor plate, and 

(d) 152 mm round bed with heater directly above. 
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The average screen diameter is found to be 

^screen x. /100 

lb- 1 ] 

d 

screen,i 

According to Levenspiel [2] , for irregular particles with no seemingly 
longer or shorter dimension, the proper mean particle size to use for fluidi- 
zation studies is 

d -±{id 

p 2 screen 

where <j> is the sphericity which can be estimated from Kunii and Levenspiel [3] 
to be $ * 0*63* Thus the average particle diameter was estimated to be 161 

jjDle 

Properties of the silicon particles were therefore assumed to be as 
listed in Table 1* 


Table 1* Properties of Silicon Particles* 


Size distribution 
Mean size 
Density* 

Bulk density 
Heat capacity* 

Thermal conductivity* 
Minimum fluidized velocity 


75 - 425 ym 
16 1 \m 
2329 kg/m 3 
1085 kg/m 3 
711 J/kg *K 
148 W/m*K 
30 mm/s 


*at 25°C, adapted from Yaws, et al [4] 


Temperature measurements were made with chrome 1 -a lumel thermocouples of 
0*254 mm diameter at locations specified in Figure 3 and were recorded using 
an Esterline Angus PD2064 data logger for steady-state tests * Transient tests 



7 
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air from 
compressor 


Fig. 3 Sketch of equipment (not to scale), location of temperature 
probes, and heat balance around the distributor plate. 
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used 127 ym diameter chrome 1-alumel thermocouples with recording of data on a 

HP 3497 data logger* 

Mr discharged from the apparatus was routed through a low pressure drop 
air washer for capture of fines which may have elutriated from the bed* 

Procedure 

The experiments reported in this paper were designed to determine the 
effects of the variation of air velocity , bed height , and power level on the 
effective configurational absorptivity and the heat transfer coefficient from 
the bed to the distributor plate* 

Two different experimental tests, steady state and transient (or pulsed) 
were used in the program* In the steady-state test, air at room temperature 
was passed upward through a bed of silicon particles at a known fixed rate to 
obtain a fluidized bed* The amount of energy input to a lamp above the bed 
was measured by a wattmeter* When the system came to thermal equilibrium, the 
inlet and outlet air temperatures, bed temperature, and distributor tempera- 
ture were measured* In the transient (or pulsed) test the bed was allowed to 
reach steady state (initially near room temperature) and then a short duration 
pulse of radiant energy was provided to raise the bed temperature approximate- 
ly 6°C above the steady-state value* The fraction of the electrical energy 
input in a pulse which was carried out by the air leaving the bed was defined 
as the effective absorptivity* 

Since temperature measurements provided the basic data for determining 
the values of h, it was important to develop a means of measuring the tempera- 
ture of the top surface of the distributor plate without disturbing the flow 
of particles, and without being influenced by the temperature of the parti- 


cles . 
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Figure 4 shows four different methods which were used to measure the top 
surface of the distributor plate® It was clear that the initial method, a, of 

fixing the thermocouple to the top of the distributor plate was not satisfac- 
tory . This thermocouple was much too responsive to the bed particle tempera- 
ture® In methods b, c, and d thermocouples were brought through the distribu- 
tor plate from below in order to eliminate the effect of the hot bed 
temperature. But in method b the temperature measurement was still affected 
by the hot particles above the distributor. For method c the temperature 
measuring point was not imbedded sufficiently in the distributor plate. 

Method d was the method finally selected. The thermocouple wires were placed 
in a shallow groove in the top surface of the distributor plate and an epoxy 
cement was used to seal the gap around the wires. 

It was observed that the distributor surface temperature fluctuated 
rapidly when air bubbles were present at the plate. As a consequence, a time- 
averaging process using 8-40 points was necessary to more accurately establish 
the mean surface temperature of the plate. 

THEORY 

Mechanism of Heat Transfer Between Distributor Surface and Fluidized Bed 

Experimentally obtained heat transfer coefficients between bed walls or 
exchanger surfaces and fluidized beds have been explained in terms of various 
mechanisms by different authors as follows: 

( 1) Steady-state conduction of heat across the gas film, which is 

scoured by solids descending along the heat exchange surface: Leva 

et al [5,6] , Dow and Jakob [7 ] , Levenspiel and Walton [8] . 





thermocouple 


Fig. 4 Different method of measure the temperature of 
the surface of the distributor plate. 
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(2) Unsteady-state thermal conduction by single particles in direct 

contact with heat exchange surface: Botterill and Williams [9], and 

Ziegler et al [ 10] • 

(3) Unsteady-state absorption of heat by fresh emulsion elements which 

are renewed intermittently by the violent disturbances in the core 
portion of the fluidized bed: Mickley et al [11,12], 

(4) Steady-state conduction through the emulsion layer which is not 
often swept away, van Heerden et al [13,14] and Wicke and Fetting 
[15]. 

It would be useful to develop criteria to suggest which mechanism con- 
trols and which type of model should be used to represent a particular situa- 
tion, Following the film-penetration theory for mass transfer, originally 
developed by Toor and Marchello [16], Yoshida et al [17] proposed a mechanism 
of heat transfer between a fluidized bed and wall surface which included both 
steady-state conduction of heat through an emulsion layer at the wall and the 
unsteady-state absorption of heat by emulsion elements. 

The mechanism of heat transfer between the distributor plate and a 
fluidized bed has not been discussed previously in the literature. It should 
be pointed out that here the direction of heat transfer is opposite to the 
direction of the air flow, rather than normal to it, as is the case in most 
heat exchangers. In comparing these mechanisms we see that these are two 
contrasting situations, as shown in Figure 5, In the first case there exists 
a film of gas at the heat transfer surface, while in the second this film is 
absent, and heat is transferred to the distributor plate primarily by impinge- 
ment of hot bed particles on the surface. 
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Fig. 5 (a) In heat transfer to bed walls or to bed internals 

a gas film is present at the surface. 

(b) In heat transfer to a distributor plate no gas 
film is present. 
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Even though the film is absent we may represent the heat transfer rate by 
a heat transfer coefficient which depends on the bed geometry , gas and parti- 
cle properties , and fluidization characteristics which in turn affect the 
collision frequency of particles with distributor plate. 

Calculation of the Effective Heat Transfer Coefficient Between the Bed and the 
Distributor Plate , h 

The heat transfer coefficient is defined so that at steady-state the 
thermal energy added to the air in passing through the distributor plate is 
equal to the energy transferred from the bed to the distributor plate by the 
action of the bed particles. Referring to Figure 3, at steady state we have, 
in symbols 



or 


m air % , air ^air , leaving ^air, entering ^ ^dis ^^bed ^dis ^ 

distributor distributor 

If the air leaving the distributor is assumed to be at the temperature of the 
top of the distributor plate then 

m . C . (T-. - T . . ) 

, air p ,air ^ dis air,m ; / 1 \ 

h A~~ [ T . T". 1 (1) 

dis ^ bed dis J 

Measuring the bed, distributor, and inlet air temperature then allows calcula- 
tion of h. 


Calculation of the Ef f ective Absorptivity of the Bed 

A certain fraction f ^ of the electrical energy provided to the lamps is 
transformed into radiant energy; and then a certain fraction f 2 of this 
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radiant energy is actually absorbed by the bed itself. The fraction f^ is a 
characteristic of the radiation source (the lamps used) , while the fraction f 2 
depends on the system geometry and the absorptivity of the bed surface. The 
product £ i @ %2 is the effective absorptivity a of the system. 

So if P is the power input to the lamps then 

^lamp to bed ^to heat air 
or 


aP 


- m . C . (T. 
air p ,air ^ 


bed 


- T ) 

air ^in' 


( 2 ) 


The factors f^ and f 2 s in particular f 2 » are the quantities which are re- 
quired; however only their product a is found directly by experiment. 


RESULTS AND DISCUSSION 

Ap vs. u Q for Three Different Bed Geometries 

Figure 6 shows the pressure drop curve for the 152 mm square bed. It is 
similar to that of an ordinary fluidized bed. 

Figure 7 shows the pressure drop-air flow characteristics for a tapered 
square bed with a 76 mm square distributor plate for three different settled 
bed heights. The velocity quoted in this figure is the superficial velocity 
based on the area of the distributor plate. It was observed that at low flow 
rates the air simply passed up through the settled solids without disturbing 
the particles. With further increase in the air flow* the surface of the bed 
suddenly erupted and there was a sudden drop in pressure. The maximum pres- 
sure drop observed is called the pressure peak. Beyond this pressure peak the 
bed bubbled and spouted * with pressure drop close to the theoretical. The 



Pressure Drop, Ap (cm H 2 0) 
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bed 


Fig. 6 Pressure drop curves for a square bed, H 


9.5 cm. 





Pressure Drop, APbed (cm H 2 0) 
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Fig. 7 Pressure drop curves for a square tapered fluidized bed 3 
distributor subtracted,, for three different settled bed 
heights . 
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pressure drop decreased progressively with further increase in velocity® No 
pressure peak was observed when the velocity was decreased® 

Figure 8 is a plot of pressure drop versus air velocity for a tapered 
square bed without distributor, or the pyramid bed, for three different 
settled bed heights® It was observed that because the gas velocity was large 
at the bottom of the bed, the particles at the bottom of the bed started 
fluidizing while those in the upper portion of the bed remained static® With 
a further increase in the flow rate of gas, layer after layer of solids became 
fluidized with the upper portion still remaining in a fixed bed state® The 
bed pressure drop thus increased continuously until the entire bed was 
spouted® Just before this stage was reached there was a sharp reduction in 
pressure drop to the theoretical® With a further increase in air velocity, 
the pressure drop decreased slightly and then increased slightly® This in- 
crease may be attributed to the energy loss by collision and friction among 
particles as well as between particles and the walls of the container® Hie 
energy loss was larger in the tapered bed than in a normal fluidized bed 
because the movement of the particles is quite different in these different 
geometries® This progression of events is consistent with other observations 
in tapered beds [18,19]® 

Both Figures 7 and 8 show that there is a minimum spouting velocity at a 
certain bed depth which becomes large as the bed height increases® 

Effect of Power Level on q 

Figure 9 shows, in the range of 300-800 W , that the effective absorp- 
tivity is independent of the power level used in the one-lamp runs® The value 
of a would be expected to change only if the lamp efficiency changed signifi- 
cantly with power level® This did not happen® 




Volumetric Flow Rate of Air, v (m 3 /hr) 


Fig. 8 Pressure drop curves for a square pyramid bed having no 
distributor s for three different settled bed heights. 
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Figure 9 also shows that a does change with bed height , hence with dis- 
tance between lamp and bed surface,, The reason is that the configuration 

factor is a strong function of the distance between lamp and bed surface* 

Effect of Power Level and Temperature Level on h 

Figure 10 shows clearly that the heat transfer coefficient in experiments 
with the resistance element located just above the bed surface and in both 
one-lamp and four-lamp experiments is independent of bed temperature for all 
bed geometries tested* 

This finding shows that the heat transfer coefficient h is a legitimate 
variable for representing the rate of heat transfer from bed to distributor 
plate* 

As a result of the above findings on power level of the lamps , the power 
level was kept at 600 W for the remainder of the one-lamp runs* 

Effect of Air Velocity and Bed Depth on h 

Figures 11, 12, and 13 show the results of experiments which examine 
these variables in beds of different geometries* 

First of all it is clearly seen, in all three bed geometries, that bed 
depth does not affect h* 

Regarding air velocity and h it is seen that in the round bed (Figure 11) 
h rises to a maximum at about 2 u m £, then decreases slowly with increase in 
gas velocity* 

On the other hand, in the square bed (Figure 12) h rises progressively 
and quite sharply with an increase in gas velocity® This is probably caused 
by the exaggerated circulation of particles in a bed of this geometry, with 
particles rising in the center of the bed, and sliding downward in the four 














24 



Fig, 13 Variation of heat transfer coefficient with air velocity for 
a square tapered bed; lamp power, P =? 600 W. 



25 


corners of the bed. This movement enhances the flow of fresh particles across 
the distributor plate, giving a higher h. An increase in air velocity will 
amplify this action causing h to increase progressively with air velocity. 

In the square tapered bed (Figure 13) h reaches a maximum at a velocity 
slightly greater than u^, drops slightly, and then gradually increases with a 
further increase in u Q . This behavior looks to be somewhere between that of 
the round and the square bed. At low air velocity in the tapered bed solid 
circulation is not severe, hence the round bed behavior dominates; however at 
high air velocity the square bed behavior with its severe solid circulation 
caused by the corners dominates. 

Effect of Air Velocity on a 

Figures 1 4 — 1 6 show the effects of air velocity on a for three different 
bed configurations. These curves clearly indicate that a increases as the 
velocity increases. This is the result of a higher average surface level of 
the bed (and therefore a larger configuration factor between the heated source 
and the bed) , more effective absorption due to the more irregular bed surface 
and the particles thrown up close to the lamp. 

It is observed that those geometries which cause more violent distur- 
bances of the bed surface result in greater rates of change of a. In terms of 
the configurations shown in Figures 2(a), (b) , and (c) ; type (c) has the 
greatest effect upon the value of a as air flow is Increased while type (a) 
has the least effect. 

Effective Absorptivity - Experiment 

Figure 17 displays the measured a values In three different kinds of 
beds, and is a crossplot of Figures 14, 15 and 16. 














ttrective ADsorptivity, a v/o) 



0 5 10 15 20 25 30 

Air Velocity, (cm/s) 


Fig. 16 Effects of air velocity on a for three different settled 

bed heights in a square pyramid bed; lamp power, P = 600 W. 




35 40 45 50 55 60 65 

Distance of lamp from settled bed surface, X (cm) 


Fig. 17 a varies with distance of lamp from settled bed surface; 
experimental values compared to the calculated values 
and to direct measurements. 
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Two observations stand out clearly: first , that bringing the lamps 

closer to the bed surface increases a, at both high and low gas velocity. 

This result is expected. Secondly * higher gas velocities give higher a 
values . This finding can be a result of two factors . 

(a) Bed expansion at higher gas velocity brings the bed surface closer 
to the lamp * 

(b) Hie severe bubbling action at higher gas velocity results in a more 
irregular surface, in which particles are thrown into the freeboard, 
much closer to the radiation source. 

We also note that at low velocity the tapered and pyramid beds have a lower a 
than the square bed. This is because their top surfaces are in tapered 
regions, and their upper cross-sectional areas are smaller than that for the 
square bed. Increasing l amplifies this effect. 

Effective Absorptivity - Prediction 

Prediction of the effective absorptivity is difficult because the geo- 
metrical effect considered in the conventional configuration factor assumes 
that radiation occurs from a diffusely radiating source. Because a major 
attempt was made to obtain a directional source rather than a diffuse one, 
this type of modeling is bound to be in error. However, the diffuse model 
should represent the lower limit for the actual situation. 

For direct radiation between diffuse source and sink Siegel and Howell 
[20] give 

f' = I [2 + 5 2 - / 5 4 + 4 S 2 J , where S = £/r 


( 3 ) 
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In the real situation the transfer of energy should be greater than that 
given by Eq . (3) * Figure 17 shows that this is so* 

Effective Absorptivity - Direct Measurement 

A series of transient tests were performed in which a pulse of radiant 
energy was directed onto a blackened copper plate (absorptivity - 1) placed at 
the surface of the settled bed* From the rate of rise of the plate tempera- 

i 

ture it was possible to get a direct measurement of the effective absorptivity 
in this bed geometry* The measured value is also shown in Figure 17* 

Note, that at close to u^-p the copper plate and the bed measurements are 
close to each other* However, at higher gas velocity where the bed expands 
(decrease in £) and has an irregular surface, with particles ejected into 
freeboard, the bed measurements give a larger a, as expected* 

SUMMARY OF EXPERIMENTAL FINDINGS AND THEIR RELEVANCE IN SCALE-UP 

L Lamp power * In all runs the effective absorptivity of the bed a and the 
heat transfer coefficient at the distributor plate h are independent of the 
lamp power used* 

2* Bed depth does not affect h* However, since bed depth changes the dis- 
tance between the bed surface and the lamp it affects the view factor, hence 
a • 

3. Bed geometry, gas velocity and q * The value of a is strongly dependent 
on the distance between radiation source and bed surface — small distance, 
large a. Also, a increases sharply with increase in gas velocity* Roughly, 
from Figure 17, we have for a: 
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u o /u mf “ 2 

u o /u mf 

l « 

48 cm 

11% 

26% 

£ - 

57*5 cm 

7.5% 

16% 


At high velocities tapered beds have higher a than nontapered beds , most 
likely because tapered beds experience more violent bubbling, more irregular 
bed surfaces, and more particles thrown into the freeboard where they are 
closer to the radiation source* 

4® Bed geometry, gas velocity and h e For all geometries tested the data 
show that at gas velocities a bit above u^ (~ 2 u^) the h values were all 
roughly about h = 260 ~ 300 W/m »K. With an increase in velocity beyond this 
point 

® h decreased slowly and progressively in round beds (see Figure 11) • 

@ h continued to rise sharply and continuously in square beds (see 

Figure 12)* 

© h decreased a bit, then increased slowly in tapered beds (see Figure 
13) • 

Thus at many multiples of \x m ^ (> 4) we would expect generally that 

h , < h , < h < h 

round round square square 

tapered tapered 

Conical and square pyramid beds have no distributor plates hence h is not 
pertinent to such systems* 
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CALCULATIONS FOR SCALE-UP 


The two conditions to be satisfied in operating the scaled-up reactor are 

that 

m the bed temperature is to be kept at about 700°C (good dense deposit 
of silicon on the particles at this temperature) * 

<§ the distributor plate should not be hotter than about 330°C other- 
wise silane will decompose while passing through the distributor 
plate and plug it a 

Base Case Calculation 

Let us start with a 1 52 mm square bed, filled to a height of 19 cm (i s 
48 cm) with particles of the size used in this study (d^ = 161 pm) , and 
fluidized with room temperature gas at a velocity of u Q = 120 mm/s (or u Q « 4 
u mf)° Let us determine the power level P and heat transfer coefficient needed 
to meet the above requirements of this process * 

1. Power input level© The required power level is found from Eq * (2), or 


aP = m C fT - T , ) 

gas p ,gas v bed gas, in' 


( 2 ) 


where 


m 


gas 


oed 


= 0 © 1 4 , from Figure 17 

= %A dis p gas = (0.120)(0.152) 2 (1.2) = 0.0033 kg/s 


C p,gas - 1004 J/kg.°C 


= 700°C 
= 20° C 


gas,iti 

Substituting all values into Eq . (2) gives 


0.14 P = (0.0033) (1004) (700 - 20) 
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from which the power level needed is 


P = 16 kW 


2. Required h . This is found from Eq. (1), or 

m C (T,. - T . 1 

, _ gas p ,gas dis gas, in' 

A dl S ^bed - T di S J 


where 


T dis 


= 350° C 
- 20° C 


gas, in 

and the other values are tabulated above . 
Replacing all values gives 


( 1 ) 


h - 


(0*0033) (1004) (350 


(0*152) z (700-350) 


20 ) ... TT/ 2 

— = 135 W/m 


Operation of the Base Case Bed 

Here we see that we need a very large power input, 16 kW, to get the bed 
up to 700°C* At higher gas flow rates the power input should increase a small 
amount, not proportionately (note that a rises sharply)* 

Since h needed = 135 W/m 2 »K while Measured = 250 ~ 300 w / m2 * K we cannot 
keep the distributor cool enough with a square or round bed* Thus we must go 
to either a tapered bed with 


A dis , 135 

A ^ 300 

surface 


0*45 


or a conical bed with no distributor, or to a multiorif ice or multicone bed 
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Scale-Up to a 450 mm Round Bed* u Q - 8 and l « 48 cm 

En ergy flux - The shape of the bed should not affect the energy flux to 

the bed* however the bed cross section and the lamp-to-bed distance should* 
From Eq« (3) the view factors for this bed size are calculated to be 

Round bed* View factor* 

152 mm 0*024 

450 mm 0 *156 


As a crude estimate let us say that a goes from 0*14 to 0*25* Thus the power 
consumption of the larger bed will be 


(16 kW)( 


450 p 0.14 
152'' 0.25 


79 kW 


Distributor plate . The heat transfer coefficient should not change much 
with scale-up* hence the same consideration 


dis 


< 0*45 


surface 

should apply to the big and the small bed* 

Increasing the gas velocity * At higher gas velocities particles are 
thrown closer to the lamps thus a is higher* But the energy flux needed is 
proportional to the mass flow rate of air* The net effect of these two oppos- 
ing effects is to require a higher energy flux to the bed, but a lower energy 
input per unit of throughput (or silicon produced)* 


Significance for Design and Final Comments 
The above calculations suggest 

1* The lamp-to-bed distance should be kept as small as possible. 
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The efficiency of utilization of radiant energy improves with scale-up „ 
The ordinary vertical sided fluidized bed resting on a distributor plate 
is not a practical design. 

The multicone or multiorifice bed is an attractive design for larger 
beds » 

Having heat flux onto the upper bed surface, rather than through the bed 
walls, and maintaining a relatively shallow bed should not lead to es- 
pecial difficulties in scale-up. 

As a possible design of this type consider the sketch below: Note that 

at the lower portion of the bed the cross sectional area increases dras- 
tically. This just happens to be the zone of rapid gas expansion due to 
reaction and due to the sharp temperature rise. 



cool product gas 

room temperature 
silane 


Taking the effective absorptivity at a = 0.2 5 the power requirement 


should be roughly (ignore the small heat of reaction and assume that 
silane and air have the same properties) 


P = m C (T - T . 1/0.25 

gas p ,gas v gas, out gas,xn' 

= (0.0033 — ^-1 ( 1004 ■-~ J - - ■) (7 00°C-350°C) 7^— 
v s ' v kg ® C ' 0 ® 2 5 

- 4.6 kW for the 152 mm bed 

2 2 

Energy flux = 4.6 kW/ (0 .15 2m) = 200 kW/m of bed surface 
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Energy Requirement = 4,6 kW/ (0 .0033 

For a feed of pure silane this translates into 

1400 kJ/kg Si produced = 0.4 kW®hr/kg Si 

and for a feed of y% silane the energy requirement is approximately 

1400/ (-j^) kj/kg Si produced - 40/y kw®hr/kg Si 

These calculations are only a crude order of magnitude approximation. At 
this stage there is little need to be more precise. 

7. Only two* nonreactive aspects of this process have been looked at in this 
preliminary research project. Other aspects of this process, such as 
© proper size distribution to use 

© stabilization of the size distribution in steady state operations 
© proper feeding and removal of product solids 

# quality of the silicon deposit in beds of different geometry 

# materials of construction for the reactive bed 

© selection of the proper radiation source, 

these should all be examined and the problems overcome, before this 
approach to the production of silicon can begin to look practical. 
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LIST OF SYMBOLS 

O 

- area of distributor plate, m 

n 

- surface area of bed, m 

- heat capacity of air, J/kg*°C 

- diameter of round bed, m 
average screen diameter, m 

- mean diameter of particles, m 

- fraction of electrical energy transformed into radiation, - 
“ diffuse configuration or view factor for direct radiation 

between two facing discs, - 

fraction of radiant energy absorbed by the surface of the 
bed, somewhat higher than f^ s - 

- heat transfer coefficient from bed to distributor, W/m *K 
settled bed height, m 

distance from lamp to surface of the settled bed, m 

- mass flow rate of air, kg/s 
heat flow rate, W 

electrical power input to the lamp, W 
pressure drop across the bed, mm 
pressure drop across a distributor, mm HpO 
temperature °C or K 
distributor plate temperature, °C 
minimum fluidization velocity, m/s 

- superficial air velocity, m/s 

- superficial air velocity based on cross section area of 
distributor plate in a tapered bed, m/s 



39 


u 86 

0 

superficial air velocity based on cross section area 
surface in a tapered bed, m/s 

of bed 

© 

V 

o 

- volumetric flow rate of air, m /hr 


x i 

- mass fraction of particles of size i in a mixture of 
cles , - 

parti- 

a 

“ effective absorptivity, - 


♦ 

- sphericity of a particle, - 


^bulk 

bulk density of particles, kg/ur 
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